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Abstract 

By combining recent HERMES data on semi-inclusive DIS (SIDIS) ^-production 
with the singlet fragmentation function , which is well determined from e + e~ data, 
we are able to extract, for the first time, the flavoured fragmentation functions , 
and without making any assumptions about favoured and unfavoured transitions. 
Whereas and are very well determined, the accuracy of is limited by 

the uncertainty in evolving from the Z° pole down to the SIDIS scale of a few 

(GeV) 2 . We discuss how the precision on could be improved. Knowledge of the 
Dq= u d s wm permit the extraction of the polarized parton densities from future polarized 
SIDIS asymmetry measurements. We study the precision that can be expected in such 
an extraction. 



1 Introduction 



Fully inclusive deep inelastic scattering (DIS) of the neutral type 



Z± + N -> l± 



+ x 



(1) 



yields information only on the combination of parton densities q(x) + q(x), in the unpolarized 
case, and on the combination of polarized parton densities Aq(x) + Aq(x) when longitudinally 
polarized leptons interact with longitudinally polarized nucleons. It is crucially from reactions 
with neutrinos and antineutrinos, in the unpolarized case, that a separate knowledge of the 
parton q(x) and antiparton q(x) densities can be inferred. This avenue is, at present, not open 
to the polarized case. 

The main approach to a separate knowledge of the Aq(x) and Aq(x) thus rests upon the 
growing activity in the field of polarized semi-inclusive deep-inelastic experiments of the type 



where h is the detected hadron. 

The cross sections (or spin asymmetries) for such reactions depend, in leading order QCD, 



fragment into hadron h. (In NLO QCD these products become convolutions.) 

It has been shown [p], @] that if systematic errors can be well enough controlled so as to allow 
a meaningful combination of data from different targets and hadrons of different charge, there 
is sufficient information to extract information on both the polarized parton densities and the 
fragmentation functions. 

In the past this has not been possible and the strategy adopted in the analysis of the 
experimental data |3|-|5| has been to assume a complete knowledge of the unpolarized densities 
q(x), q(x) and of the fragmentation functions D^(z), D^(z). With these, in ||] an auxiliary 
function was constructed, the flavour (g = u,d, s,u,d,s) purity P^ N (x,z) |J for each hadron 
h and for each target nucleon N. Given the purities, the polarized data can then, in principle, 
be used to directly extract the polarized densities Aq(x), Aq(x). 



(2) 



upon products of parton densities and fragmentation functions 
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However, the situation has now changed because recently the HERMES group has for the 
first time published unpolarized charge separated data for ^ production on a proton target 
0. The main aim of our paper is to demonstrate that this data, taken in conjunction with the 
information on the flavour singlet combination of FFs which can be fairly reliably obtained 
from the data on e + e~ — > ir X at the Z° peak, allows a first direct determination of the FFs 
, D1 + and -DJ + . We note the caveat that the data in |7] covering 0.2 < z < 0.9 exhibit large 
0(40%) isospin violations at large z > 0.7. It is most unlikely that such a large breaking of 
isospin invariance can be a genuine effect in the current fragmentation region of semi-inclusive 
DIS (SIDIS), and it would seem unnatural to incorporate it into an FF formalism. A possible 
explanation for the effect is given in 0. Our analysis will therefore be relevant mainly for 
intermediate 0.2 < z < 0.7. 

It turns out that the least well determined FF is -DJ + , since it is most dependent on the 
evolution downwards of (z, Q 2 ) all the way from Q 2 « m 2 z0 to Q 2 « few (GeV) 2 , and this 
involves mixing with the gluon FF D^(z,Q 2 ). We try to assess the sort of accuracy required 
in future e + e~ measurements in order to achieve an accuracy of 20-30% on (z). 

In order to study the accuracy of the polarized parton densities obtained in the past via 
the use of purity [|J we proceed as follows. Firstly, we use our FFs to calculate the central 
value and errors of the integrated purity function J dz P£/ p (x,z) and j dz P^ p (x,z) for a 
proton target. Then, because at present the polarized SIDIS data does not exist separately for 
7r + and 7r~, we take the central values of the published polarized parton densities of Leader, 
Sidorov and Stamenov || obtained purely from DIS data, and using the central values of the 
purity we generate fake "data" on polarized SIDIS asymmetries AA™ + and AA^ , and also on 
the polarized DIS asymmetries. We then go through a similar procedure as adopted by the 
HERMES group to obtain from this "data" the polarized parton densities, with this difference, 
that we allow for the uncertainty in the value of the purity function. In this way we obtain an 
indication of the uncertainty in the polarized parton densities inherent in the purity approach. 

In this paper we work to leading order (LO) in QCD as the purity concept only makes sense 
in LO and because in LO we can deal with simple algebraic equations which are physically most 



2 



transparent to interpret and have a well-defined error propagation. Of course, in the long run, a 
more complete NLO analysis will be required. Standard experimental techniques || are based 
on an ad hoc combination of LO (polarized) parton distributions with e.g. LUND-type Monte 
Carlo fragmentation functions. Such an effective approach cannot be extended to NLO without 
a highly non-trivial definition of the long- and short distance pieces in the MC environment 
which is - to our knowledge - lacking at present. It will, therefore, be vital to bring the 
measurements in touch with well-defined factorizecf] QCD approaches |J combining universal 
parton distribution functions (PDFs) with universal FFs because, otherwise, the extracted 
PDFs will not have any physical relevance. The most important experimental information 
will be on scheme- and mo del- independent data for cross sections and not on the extracted 
(unobservable) PDFs and FFs. 



2 Extraction of Fragmentation Functions 



2.1 Formalism 

For a leading order treatment we follow the notation of |2| and remove some kinematical factors 
by introducing for the DIS and SIDIS cross sections on a proton target: 



-DIS 



-DIS 



a" 



Act 



x(P + l) 2 

4ira 2 
x{P + l) 2 

x{P + l) 2 
4na 2 

x{P + l) 2 
4na 2 



2y 2 



1 + 
y 



2„DIS 



d 2 a 



y 
2y 2 



y) 2 J dxdy 
d 2 a°L s 



d 2 a5 IS 



dx dy dx dy 



1 + 

y 

2-2/ 



d?o h 



y) 2 J dxdydz 
d 3 a^ d 3 



(3) 
(4) 
(5) 
(6) 



dx dy dz dx dy dz 

Here, P M and Z M are the nucleon and lepton four momenta, and o\ v refers to a lepton of helicity 
A and a nucleon of helicity v. The variables x, y, z are the usual DIS kinematic variables. Then 



lr The factorization of x— and z— dependence at LO is an artefact from the one particle phase space (delta- 
function) of LO-SIDIS and is of no fundamental importance as opposed to the mass factorization of NLO cross 
sections. 
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one has the very simple LO results: 

a DIS (x,Q 2 ) = 5>^(*,Q 2 ) (7) 

Aa DIS (x,Q 2 ) = ]Te 2 Ag 4 (*,Q 2 ) (8) 

9:5 

Aa h (x,z,Q 2 ) = ^e 2 Ag 4 (x,Q 2 )Z^(z,Q 2 ) (9) 

9:5 

a\x,z,Q 2 ) = ^e 2 % (x,g 2 )^(^g 2 ), (10) 

9:5 

Note that the inclusion of a factor (1 + R)/(l + 7 2 ) in (|10D (see e.g. Eq. (5) of [f|) is 
not justified theoretically. The correct handling of the longitudinal cross-section is a more 
complicated NLO effect in SIDIS (see Eqs. (56) - (60) of 0). Here, as mentioned, we work to 
LO only. Specializing now to 7r ± production we introduce the measured observablesQ 

Rf(x,z,Q 2 ) = (11) 
v 



ft. 

a Dls 
v 



:i2i 



Using charge conjugation and isospin invariance we require only 3 independent FFs: 

Df(z,Q 2 ), Df(z,Q 2 ), Df(z,Q 2 ) (13) 
The remaining ones are then: 

(14) 
(15) 
(16) 

Thus 

= ^Js { \ (uDZ + + uDt) + I (dDf + dDf + sDf + wf) } 

= -L^I(Au + d)D: + + (Au + d)Df + (s + s)Df} (17) 





= Df 


= Df 


= Dl 


Df 


= Df 


= Df 


= Dl 


Df 


= Df 


= Df 


= d: 



2 As we are considering positive and negative charges seperately we are not adopting the convention = 
h+ + h~ ilMl. 
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Similarly 

R T = ^dis { ( m + d ) D u + + ( 4m + d)Df + (s + s) Df } (1 
Assuming a good knowledge of the unpolarized parton densities we can immediately obtain 



9 [Rp - Rp )a 

In order to obtain + DJ and £)J we require one further piece of experimental in- 
formation. We shall argue that it can be obtained from the data on e + e~ — > ir X at the Z° 
peak. 

2.2 Use of the e + e data 

For some time it was believed that the fragmentation functions obtained by Binnewies et al. 



10], from a detailed analysis of the e + e data over a wide range of energies, were reasonably 



well determined. However, recent analyses [[Tl| -|T3|j have shown that equally good fits to e + e~ 
data can be achieved with FFs of a given flavour which differ widely from each other. The e + e~ 
data do not, therefore, constrain the FFs of a given flavour very well, and, in retrospect, this 
is really not surprising. 

However, by a piece of good fortune, the e + e~ data at the Z° peak directly measure a linear 
combination of FFs which is very close to the SU(3)f flavour singlet combination, i.e. in 

D£2T = E (D n q ++W ~ + Df +n ~) e 2 q (s) (20) 

q=u,d,s 

the squared electroweak couplings e 2 (s) from SU(2) x U(l) gauge symmetry (given e.g. in the 
appendix of [jTTH ) are flavour- independent to within ~ 25 % at yfs = Mz as opposed to a 
relative factor of 4 for the electromagnetic couplings of up- and down-type quarks at lower cms 
energies. The exact singlet below in ( p2|) would correspond to a measurement at an e + e _ cms 
energy of = 78.4 GeV or y/s = 113.1 GeV where it happens that e^(s) = e^(s) = e%(s). 
Accordingly (|20| ) is approximately proportional to the singlet combination 

D£ = (Df + Df + Df+Df+Df+Df) (21) 
= 2(Df +Df +Df) (22) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

z z 



Figure 1: The ratio of singlet fragmentation functions D% found in |TT| (K) and [0 (KKP), 
at Q 2 = 100 GeV 2 (left) and at a typical SIDIS value Q 2 = 2 GeV 2 (right). 



where we have used charge conjugation and eqs.([L4]) - flIE|) in the last step. Using isospin and 
charge conjugation invariance Z?£ +D? = +DJ ± and approximating e 2 1 (s)/e^(s)| s=A/2 = 
G-u(s) /e 2 (s)\ s=M 2 - 3/4 we can write the singlet combination 

D£ = ^&£+T - \ (Of + Df) (23) 



ij me as y 

where we have introduced a convenient change in normalization 

= D^tr I e 2 {s) (24) 
The extreme limits < (D^ + + D? + ) < (D^ + + D? + ) then correspond to 

\ < D£ < i Dtt:~ (25) 

i.e. to only a ~ 5% uncertainty for . 

Not surprisingly, therefore, the singlet FFs in the analyses fil"|-|i"3"[ agree with each other to 
better than 5% for 0.2 < z < 0.7 as seen in Fig. |l|. We may thus take as a known quantity 
D^ + (z, Q 2 = m 2 z0 ) and from Fig. [I] we observe a stable evolution down to Q 2 = 100 GeV 2 . 

But we require this quantity at a scale of a few (GeV) 2 and it thus has to be evolved down 
through a large range of Q 2 , and in this evolution mixes with the poorly known gluon FF D 1 ^ . 
(Of course we cannot carry out the evolution of -D mea s itself since it is a combination of singlet 
and non-singlet pieces and we do not know the values of these separately.) The FF analyses in 



10|- |i~3[l cover data down to \fs ~ 30 GeV and from Fig. [I] we judge this fixes a stable singlet FF 



down to y/s ~ 10 GeV. Below, however, the evolution uncertainties set in and from the right 
of Fig. [l] we quantify this uncertainty conservatively to be a ~ 20 % effect uniformly in z. We 
convinced ourselves this is indeed a typical order of magnitude by comparing the several sets of 
LO and NLO FFs for tt+ + n~ , K + + K~, h + + h~ in and not only the two sets plotted 



in Fig. |l|. Clearly, a low scale measurement of the singlet FF or a resolution of the evolution 
ambiguities through a determination of the gluon FF would be highly desirable information. 

Subject therefore to possible errors due to the evolution, we have available the additional 
experimental data that we require, and we then obtain: 

9(R; + + Rf)a° IS -2sDl + 
4(ii + u — s) + d+ a 

and 

-18 (Rf + Rf~) d° IS + [4(u + u) + d + d] Dl 

jy* = \— !— = f27) 

2[4{u + u-s) + d + d] v ; 

We note that the singlet FF plays no role in eq. ([19]) and that its weight increases in going from 

(|26|), where it is multiplied by the suppressed strange PDF, to (|2T|), where it is multiplied by 

the unsuppressed u and d PDFs. Correspondingly, we must expect an increasing importance 

of the propagation of the ~ 20% error of into these equations. 

The LH sides of eqs. (p!9|),(p6|) and are functions of z and Q 2 , whereas the RH sides 

are, in principle, functions of x, z and Q 2 . Only in the LO approximation does the variable 

x become a passive variable i.e. there is no dependence on it. Strictly one should test for 

this lack of x-dependence, as a measure of the reliability of the LO treatment. However, in 

this paper, in order to improve statistics, we shall take it for granted that the LO treatment is 

adequate. 

2.3 Combined analysis of HERMES and e + e data 



The formalism given in (|14])-([16|), (|i9|), (f22j), (|26|) and (|27|) presupposes the availability of data 



at fixed x and y. In fact the available HERMES data is integrated over the kinematic range 
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Q 2 > 1 (GeV/c) 2 , W 2 > 10 GeV 2 , y < 0.85. Handling integrated data slightly complicates the 
formalism since the ^-dependent factors in the numerator and denominator of (JTTJ) no longer 



cancel out to give the simpler result (0). We have done the analysis using the data integrated 
over the kinematic range of the experiment and have checked that using the simpler formalism 
with gj 

x =< x >= 0.082 (28) 
Q 2 =< Q 2 >= 2.5 (GeV/c) 2 (29) 
W 2 =< W 2 >= 28.6 (GeV) 2 (30) 

makes no discernible difference to the results for the FFs. 

The stability of our results for the central values of D^ + , is studied in Fig. The 
NLO determination of the singlet combination due to Kretzer |l]J was utilized. To test the 



stability of our results we have used two different sets of unpolarized parton densities. We found 
the effect of employing, respectively, the NLO MRST |TJ[ or the NLO GRV [[15] unpolarized 
parton densities in eqs. (0), flUf), (|19"D , (ffif) and ( p7f) leads to a negligible ~ 5% effect. This is 
because the unpolarized densities are very well constrained in the region of interest. We have 
checked that use of the LO GRV densities also have almost no noticeable effect. 

The FFs D% and Df are quite well constrained by the SIDIS data, but, as expected, 
is sensitive to the singlet combination of FFs determined from e + e~ data. This can be clearly 
seen in Fig. ^| where we compare the results for each + using NLO versions of as obtained 



by Kretzer [1 1 and Kniehl, Kramer and Potter [12] from the e + e~ data. As mentioned in the 
Introduction, is very well determined at the Z° peak, but the mixing, under evolution, 
with the gluon FF induces an uncertainty of about 10 — 20% at Q 2 = 2.5 (GeV/c) 2 . As 
can be seen in Fig. |^, may even turn unphysically negative at large z > 0.7 where our 
analysis is not supposed to be reliable, anyway, as mentioned in the Introduction. We note 
Fig. |^ shows the typical effect expected from the uncertainty of and that a similar picture 
emerges if we switch between the LO and NLO parametrizations of (jllj: The stability of 
and DJ + is remarkable. -DJ + , on the other hand, changes significantly. If a measurement 
could fix (Q 2 ~ 2.5 GeV 2 ) to within ~ 5% we would have a handle on as well at the 



h- 
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Figure 2: The valence- type (D^ + ), sea-type (-DJ + ) and strange- type (Dg + ) FF into charged 
pions extracted under the assumption of isospin-invariance from HERMES SIDIS measurements 
supplemented by the singlet FF of [JTTJ] (solid) and |12| (dashed), respectively. Details are given 
in the text. Switching between the LO <-» NLO parametrizations of ]TT],|T^[ leads to similar 
variations in the FFs whereas the variation from employing different unpolarized PDFs is 
negligible. 



~ 20 - 30% level. 

In summary we see that and are remarkably well constrained by the SIDIS data. 
_DJ + however, is undetermined within a factor of about 2. In Fig. |^ we show the final results 
for our FFs and our estimate of their errors. These include Gaussian error propagation of the 
experimental errors in H combined with a 20 % error of D^(Q 2 = 2.5 GeV 2 ). The FFs can 
be described analytically 3 at (Q 2 ) = 2.5GeV/c 2 by 

Dl + = 0.689 z- img (1 - z) 1 - 241 (31) 

3 Employing D™ + — N z a (1 — z)^ q ansatze with flavour independent N,a at (Q 2 ) — 2.5GeV/c 2 slightly 
worsens the quality of the parametrization. 
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Figure 3: The extracted fragmentation functions with errors which combine the experimen- 
tal errors from Q with a typical 20% uncertainty arising from the evolution of the singlet 
fragmentation function. 
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Figure 4: A comparison of the extracted FFs , (solid lines) at (Q 2 ) = 2.5 GeV 2 as in 
Fig. |l]to the LO-parametrization in |Tl| (dashed) where (1 — z)D^ + (z) = (z) at the input 



Ql 



Df = 0.217 z" 1 - 805 (1 - z) 2mr 



D1 = 0.164 z 



-1.927 



(32) 
(33) 



where _DJ + has to be taken with care as Fig. |] shows that basically no value within < _DJ + < 
can be excluded at present. 
Finally, as a matter of interest, we compare in Fig. f| the D^ uds obtained in this paper 
with the LO obtained by Kretzer |llj purely from an analysis of the e + e~ data. In the 
latter the flavour separation is not fixed by the data and is somewhat ad hoc and assumed 
D% + = Dg + and > by imposing (1 — z) = DJ + at the input scale. Surprisingly, 
Kretzer's are not very different from the obtained from our analysis of the SIDIS 
data! 
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3 Implications for the polarized parton densities 

As mentioned in the Introduction, the absence of neutrino data for polarized DIS means that 
the extraction of the individual Aq(x) is impossible. Only the combination Aq(x) + Aq(x) can 
be found and the flavour separation of these relies heavily on the evolution in Q 2 and is thus 
unreliable, given the small range of Q 2 available in polarized DIS experiments. Thus polarized 
SIDIS has a vital role to play in this matter. 

At present, however, there are no published data for polarized ^-production, though there 
does exist data for undifferentiated polarized h^- production, which have been used by the 
HERMES group to extract information on the polarized parton densities via what is known as 
the purity method. 

We believe that in this approach the errors on the polarized parton densities are somewhat 
underestimated and we shall use our FFs to study this question. 

The flavour q purity function for protons || used by the HERMES group f||] is defined by 

Jo.2 D q( Z ) dZ 

E^^q'ix) Sl 2 D h q ,{z)dz 

where again, we utilize the MRST parton densities and take Q 2 =< Q 2 >. 
Defining now the SIDIS spin asymmetry 

, f„„ dz Aa^ix, z) 
< AA h p (x) >= ?—L2. (35) 

J 02 dzcr$(x,z) 



we have in LO, 

' Aq(x) 



< AA h p (x) >= (" 



(36) 



q(x) 

Similarly for the DIS spin asymmetry we can define 
& 2 q{x} 

Pq/p i X ) — 2 TT~\ ffl) 



4 Note that the graphs shown in the HERMES publications |16|] and labelled "purity" are actually plots of 
an "effective purity" incorporating various experimental cuts. 
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and then, in LO, we have, with Q 2 =< Q 2 >: 

Similar expressions for < AA^{x) > and Aa^ IS (x) / a^ 1 (x) can be obtained in an obvious way. 

At each value of x there are in principle 6 pieces of data (h = 7r ± for p, h = for n, and 

DIS for p, n), so that there is enough information to solve for the 6 quark polarized densities 

Aq(x) I q(x), for q — u,u, d, d, s, s. In the published analyses of the latter data || the HERMES 

group has preferred to model the polarized sea with assumptions such as 

Am Ad As As Aq s . . 

— = = — = — = ^i± 39 

u d s s q s 

or 

Au = Ad = As = As = Aq s (40) 

and then to obtain the 3 independent polarized densities by making a best fit to the 6 -pieces 
of data at each x. 

The problem with this approach is that the purity functions were constructed using LUND 
model information on the FFs. We think [|2[ this is an un-reliable procedure since a combina- 
tion of (polarized) PDFs and LUND-type of FFs is at present lacking a rigorous theoretical 
framework as opposed to our combination of universal (polarized) PDFs with universal FFs in 
line with the factorization theorems of QCD []17|| . 

We argue that the above approach much underestimates the uncertainty on the polarized 
parton densities. To illustrate this we construct purity functions and their errors for pion 
production, using the fragmentation functions determined by us and the unpolarized MRST 
parton densities. The formulae are exactly as in fl34|) , (|35l) and (|36|) with h replaced by 7r + . We 
show our calculated purities for protons with errors in Figs, ^j, |6|. The size of the errors are 
as to be anticipated from the errors on the DZL u ds m Fig. [3| and the definition of the purities 
in Eq. fl34|) . In the absence of separate 7r ± SIDIS spin asymmetry data, we take the central 
values of the polarized parton densities, as derived by Leader, Sidorov and Stamenov M from 
purely DIS data and by feeding these Aq into eqs. (^) and ( |3"%D and into the analogous one 
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Figure 5: Purity functions Pq/ p {% 
from the evolution of the flavour- singlet D^. 



10 10 

extracted from HERMES data 



with 20% uncertainty 



for AAp , in which we utilize the central values of the purity functions, we generate a set of 
fake "data" for AAf , AAf and Aa° IS /a° Is . 

Having now this set of fake "data" we forget where it came from and use it to solve for the 
polarized parton densities, mimicking the approach used by the HERMES group. Thus we take 

Au Ad As Aq s . . 

— = -=- = — = 41 

u a s q s 

and solve (p8[) , (|35|) (for h = 7r + , 7r _ ) for Au/u, Ad/d and Aq s /q s . In this analysis we treat the 
"data" as perfectly known, but include realistic errors on the purities, arising from the errors 
on our FFs. In this way we illustrate the uncertainty on the polarized parton densities arising 
solely from the uncertainties on the purity functions. 

The results are shown in Fig. [7[ It is seen that whereas Au/u is largely insensitive to the 
uncertainty on the purity, both Ad/d and Aq s /q s inherit significant errors from this uncertainty. 
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Figure 6: Purity functions Pq/ p (x 
from the evolution of the flavour- singlet D^. 



10 10 

extracted from HERMES data |7]] with 20% uncertainty 



Bearing in mind that the errors shown in Fig. [7] arise solely from the uncetainty on the 
purities, one learns from this study that it is misleading to treat the purities as absolutely 
known quantities. It would be far more meaningful to follow the strategy suggested in || and 
use the SIDIS data to obtain both the FFs and the polarized parton densities. The purity is 
an unnecessary element and in any case loses its usefulness in NLO. 



4 Conclusions 

We have shown that a judicious combination of the HERMES SIDIS data on ^ production 
and certain aspects of the data on e + e~ — > if^X allows the extraction, for the first time, of the 
flavour separated fragmentation functions D* + , DJ + and D* + . 

The key element in this approach is the avoidance of any ad hoc model-dependent flavour 
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Figure 7: Polarized parton densities extracted from the perfect "data" as explained in the text. 
The HERMES simplifying assumption Au/u = Ad/d = As/s = Aq s /q s has been used. The 
uncertainties arise solely from the realistic experimental errors on our fragmentation functions. 
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separation in the e + e data, by noting that at the Z° peak what is well determined is essentially 
the light flavour singlet combination of fragmentation functions _D£ + = 2 ( + + ) , 



which is almost identical in all analyses of the e + e~ data. The negative aspect of this approach 
is the need to evolve -D£ + down from the Z° region to the SIDIS region of a few (GeV) 2 , which 
involves mixing with the poorly known gluon fragmentation function. 

In fact the HERMES data is extremely accurate, so that almost all the uncertainty in our 
determination of , -DJ + and D J + arises from the uncertainty in the evolution of _D£ + . As it 
turns out, this has little effect on and , which are very well determined, but has 
relatively large errors. 

We have also examined the question of the precision with which the polarized parton den- 
sities could be extracted from future polarized SIDIS pion production data. Here we have 
assumed perfect 'data', then followed the HERMES purity method to obtain the polarized 
parton densities, and thereby displayed the uncertainties generated solely by the errors on the 
purity functions. The significance of this study is that in the earlier analyses [Q, |j] the purity 
functions are taken as almost perfectly known with essentially no errors. As expected we have 
found that Ad and Aq s are significantly affected by the uncertainties in the purity functions. 
This suggests that in the published polarized parton densities extracted from polarized SIDIS 
/^-production data |3|-[5| , the uncertainties given are missing an inherent error arising from the 
fragmentation uncertainties as quantified in this paper. 
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